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Asymptotic  Probabilities  in  a 
Critical  Age-Dependent  Branching  Process 

by  Howard  J.  Weiner 


ABSTRACT 

Let  Z(t)  denote  the  number  of  cells  alive  at  time  t in  a critical 
Bellman-Harris  age-dependent  branching  process,  that  is,  where  the  mean 
number  of  offspring  per  parent  is  one.  A comparison  method  is  used  to  show  for 
k > 1,  and  a high-order  moment  condition  on  G(t),  where  G(t)  is  the  cell  lifetime 
distribution,  that  lim  t P[Z(t)=k]  = »k  > 0,  where  {ak}  are  constants. 

t -*  09 

The  method  is  also  applied  to  the  total  progeny  in  the  critical 
process . 
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Asymptotic  Probnbi  1 i t i c.b  In  a 
Critical  Age-Dependent  Branching  Process 


by  Howard  J.  Weiner 


1.  Introduction. 

Let  Z(t)  denote  the  number  of  cells  alive  at  time  t in  a critical 
Bellman-Harris  [2]  age-dependent  branching  process.  That  is,  at  time 
t =•  0,  the  process  starts  with  one  new  cell  with  lifetime  distribution 
G(t),  G(Of)  = 0,  non-lattice,  and 


(1.1)  G'(t)  = g(t), 
the  density  exists,  and  for  some  6 > 0, 

r*  4+6 

(1.2)  t g(t)dt  < «*. 

J0 


At  the  end  of  its  life,  the  cell  disappears  and  is  immediately  replaced  by 


k new  cells  with  probability  p^  > 0,  with 


(1.3)  2 Pk  - 1 

k=0 


and  criticality. 


(1.4)  2 kp.  - 1. 

k-1 

Each  new  cell  proceeds  identically  as  the  parent  cell  and  independently 
of  all  other  cells  and  the  state  of  the  system.  Assume  that  for  all  integers 

" k 

(1.5)  T.  up  < ». 

n-1 


It  la  the  purpose  of  thla  note  to  show  that  for  k>  1, 

(1.6)  lim  t2P[Z(t)=k]  - *.  > 0, 

t 

where  {a^}  are  constants. 
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The  method  Is  also  applied  to  total  progeny  In  the  critical 
process . 

2.  Integral  Equation  and  Result. 

Define  the  offspring  generating  function,  for  0 < a < l, 

00 

(2.1)  Ms)  3 £ Pk*  • 

k=*0 

Define  the  generating  function  for  Z(t)  as,  for  0 < s < 1,  *11  t > 0, 

(2.2)  F(s,t)  3 E(sZ(t)). 

Then  [2]  F(s,t)  satisfies 

(2.3)  F(s,t)  - 1 - G(t)  + > h(F(s,t-u))dG(u), 

- 0 

with  F(s,0)  =■  *• 

* 

Theorem.  Let  Z(t)  denote  the  number  of  cells  »t  time  t in  a critical 
Bellman-Harris  branching  process  satisfying  (1.1)  - (1.5).  Assume  also 
that  the  derivatives  of  h satisfy 

(2.4)  0 < h"(l)  < -,  0 < *»'"  (1)  < •• 

Then  for  k > 1 

(2.5)  11*  t2P[Z(t)-k]  - *k>  0 
t -»  • 

where  (a^J  are  constants. 

Proof.  Consider  first  the  case  k ■ 1,  and  write  for  simplicity 
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(2.6)  P(t)  a p[Z(t)=l). 

From  the  representation  (2.3),  derivatives  of  F(s,t)  with  respect  to 
s exist,  and  note  that 

(2.7)  — = p(t)* 

08 

From  (2.7),  (2.3)  one  obtains  the  relationship 

(2.8)  P(t)  =■  1 - C(t)  + J h' (l-Q(t-u))P(t-u)dG(u) 

0 

and  P(0)  =*  1, 


where 


(2.9) 


From  [2] , 


(2.10) 


where 


(2.11) 


Q(t)  = P[Z(t)  > 0J. 


lim  tQ(t)  - 2^(h"(l))_1  , 

r — • oo 


r 

I tdG(t). 
*■  0 


Assume  for  simplicity  in  all  the  following  that 


(2.12) 


H s i. 
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The  proof  will  proceed  by  a number  of  claims. 

Claim  I.  Let  R(t)  be  continuous,  R(0)  = 0,  and  satisfy,  for  t > 0, 

that 

(2.13)  R(t)  > f h'(l-Q(t-u))R(t-u)dG(u). 

«)  Jo 

Then  for  all  t > 0, 


(2.14)  R(t)  > 0. 

«) 

Proof.  Assume  the  upper  inequality  of  (2.13). 

Note  first  that  for  all  t > 0, 

(2.15)  R(t)  * 0. 

This  follows  from  (2.13)  by  assuming  there  is  a t^  such  that  R(tQ)  * 0, 
and  a contradiction  is  clear. 

If  the  upper  Inequality  in  (2.13)  Is  false,  (2.15)  requires  that 
R(t)  < 0 for  t > 0.  Then  on  an  arbitrary  interval  [0,t],  assume  that  R(t) 
assumes  its  minimum  at  t^  < t.  Since  0 < h’(x)  < 1,  setting  t * t^  in  the 
upper  inequality  (2.13)  yields  a contradiction.  The  lower  inequality  of 
(2.13)  Implies  the  lower  inequality  of  (2.14)  similarly,  completing  the  claim. 


Claim  II. 
(2.16) 


For  some  0 < c < 1,  all  t > 0,  define 


t < 1 
t > 1 


» 


t 


(2.17) 


L(t) 


t < 1 
t > 1. 
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Then  for  all  t sufficiently  large. 


(2.18)  K(t)  > l - G(t)  + f h'(l-Q(t-u))K(t-u)dG(u) 

0 


and 


(2.19)  L(t)  < 1 - G(t)  + I h'(l-Q(t-u))L(t-u)dG(u). 


Proof.  To  show  (2.18),  note  that  the  right  side  yields  the  Inequality 


(2.20)  I h' (l-Q(t-u))K(t-u)dG(u)  < G(t)  - G(t/2) 

J0 

pt  / 2 

+ h'(l-Q(t-u))K(t-u)dG(u). 

J0 


A Taylor  expansion  of  the  rhs  of  (2.20)  yields. 


rt/2 

(2.21)  J h' (l-Q(t-u))K(t-u)dG(u) 


pt/2 

- I (l-Q(t-u)h"(l)  + o(Q(t-u)))K(t-u)dG(u) . 
1 0 


A Newton  binomial  expansion  applied  to  the  rhs  Integrand  of  (2.21) 


using  (2.10),  Theorem  4,  p,  406  of  [3],  (2.16),  and  that  p.  s 1 yields  that 
the  rhs  of  (2.21)  equals 


6 


<t/2 


<*•*>  i 


(t-u) 


1 rt/2  dG(u) 
t2'C  0 ..  ux-  - t 


(1  -V 


2 . , 1 y 

2-e  ' 3-e  + o(.  3-e) 

ux  t t 


,t/2 


b L (1  + - ~b  + °<-b> 


+ s&si . 

1 o 


t2_€  t3-6 


2 . , 1 „ 
7^  + 


2-e 


for  alt  t sufficiently  large.  In  view  of  (2.20),  (1.1),  this  suffices 
for  (2.18),  and  (2.19)  follows  similarly. 


Claim  III.  For  0 < e < 1,  and  all  t sufficiently  large, 

(2.23)  -rrr  < p(t)  C . 

t t 

Proof.  This  follows  from  Claims  I and  II,  where  R *>  K - P,  for  example. 
Claim  IV.  Let  R(t)  satisfy 

(2.24)  R(t)  - 1 - G(t)  + f h' (l-Q(t-u) )R(t-u)dG  (u) , 

J0 

with  R(0)  ■ 1,  where  Q(t)  Is  the  probability  of  non-extlnctlon  for  the  process 
with  lifetime  G(t),  and 


1 - Go(t)  - e_t. 


(2.25) 
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Then  as  t -*  00 , 

(2.26)  R(t)~~ 

t 

for  some  constant  c > 0. 

Proof.  Write  (2.24)  as 

»t 

(2.27)  R(t)  - 1 - G(t)  + e | h’ (l-Q(u))R(u)eudu. 

1 0 

Differentiating,  with  respect  to  t, 

(2.28)  R'(t)  = -g(t)  + (l-G(t)-R(t))  + h'(l-Q(t))R(t). 

Then,  expanding  h in  a Taylor  series,  one  obtains 

(2.29)  R'(t)  +|  R(t)  - f(t) 

where 

(2.30)  f(t)  = 1 - G(t)  - g(t)  + -yV(t)R(t)  + (Q(t)  - h»(l)^)R(t>» 

with  1 - Q(t)  < * < 1.  In  view  of  (1.1),  (1.2),  (1.5),  It  follows  from 
Theorem  4 of  [3],  p.  406,  that 

(2.31)  (Q(t)  - - 0(-i2|-£). 
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An  analysis  similar  to  that  in  Claims  I - III  yields  that  for 
0 < e < 1,  for  t sufficiently  large, 


(2.32) 


-is<R(t)  <-^j 


Then  (1.1),  (2.31),  (2.32)  used  in  (2.30)  imply  that 


(2.33) 


t2f(t) 


is  integrable  with  respect  to  Lebesgue  measure  on  the  positive  line. 
To  solve  (2.29)  for  t large,  let 


(2.34) 


Q(t)  = tZR(t). 


Then  (2.29)  becomes 


(2.35) 


Q'(t)  - t*f(t), 


For  fixed  t >0,  one  then  obtains 
o 


(2.36)  Q(t)  - Q(to)  - J S2f(§)dS 


(2.37) 


R(t)  ~ 

t 


for  t sufficiently  large,  and  c > 0,  completing  this  claim. 
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Claim  V.  let 


„t 


(2.38)  T(t)  = 1 - G(t)  + j h'(l-Q(t-u))R(t-u)dG(u) 


0 


Then,  as  t -•  00 , 


(2.39) 


t2|T(t)-R(t)|  - 0. 


Proof. 


t/2 

(2.40)  |T(t)-R(t)|  < |J  h'(l-Q(t-u))R(t-u)(dG(u)-dGo(u) 


+ G(t)  - G(t/2)  + Go(t)  - Go(t/2). 


A Taylor  series  and  Newton  binomial  expansion  on  the  right  side  integrand 
of  (2.40)  yields 

(2.41)  j 1 l-Q(t-u)h"(l)  + Q2(t-u)h"'  (k(t-u))"T + o ( -~)  j'dG(u)-dC  (u> 

J0  L JL(t-u)  (t-u)ZJv  ° 


where 


(2.42) 


1 - Q(t-u)  < k(t-u)  < 1, 


and  expression  (2.41)  equals 
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(2.43) 


(t-u) 


+ L-^)J  (dG(u)-dGo<u)) 

(t-u) 


|(1  + p + o{±)]  + o(-i)]  (dG(u)-dGo(u)) 


-4 

= o(t  *), 


for  t sufficiently  large,  using  that 


H = J udG(u)  =*  J udGQ(u)  = 1. 

This  completes  Claim  V. 

The  proof  of  the  theorem  for  P(t)  may  now  be  completed. 
Define  the  iterative  scheme,  for  n > 0 


(2.44)  P(n+1)(t)  = 1 - G(t)  +|  h'(l-Q(t-u))P(n)(t-u)dG(u), 


and  let 


(2.45)  P(0)(t)  5 R(t)* 

Then  from  (2.38),  it  follows  that 

(2.46)  P(l)(t)  " T(t)> 

Claims  IV  and  V imply  that,  for  t large, 

(2.47)  • P(Q)(t)~-^  and  P(i)(t)~-^. 


j 


11 


An  Induction  argument  along  the  lines  of  Claim  V establishes  that, 
for  large  t, 


(2.48) 


Note  that  for  n > 0, 

pt 

(2.49)  P(t)  - P(nfl)(t)  -J  h'(l-Q(t-u))(P(t-u)-P(n)(t-u))dG(u). 


Denote 

(2.50)  An(t)  = P(t)  - P(n)(t), 


then  from  (2.49) , 

(2.51)  lAtvfl<t)|  < J lA„(t-u)|dG(u)  < lAol^Ct). 

where  * denotes  convolution,  and  Gn(t)  is  the  nth  convolution  of  G with 
itself . 

Let  {x^}  denote  I.I.D.  random  variables,  with  distribution  function  G(t). 

Set 


(2.52) 

Then 


n 

E X. . 
i-1 


(2.53)  Gn(t)  3 P[Sn  < t]  - P[Sn-n  < t-n] , 


w 


r 
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and  if  one  takes  n > t,  by  (1.2)  and  Chebyshev's  inequality, 


Var  S 


(2.54) 


G„(t)  < 


n n 


2 2 * 
(n-t j (n-t)4 


Also, 


t/2 


(2 


.55)  |Aq|  * Gn(t)  < G(t)  - G(t/2)  + J |p(t-u)-R(t-u)|dGn(u) 


and  (1.2),  Claims  III,  Iv,  and  (2.54)  yields  that 


(2.56) 


i4oi  • 


e*  ‘(n-t)' 


where  K is  a constant. 

For  fixed  large  t,  let 


(2.57) 


n > [t1+€]. 


This  suffices  for  the  proof  of  the  theorem  for  P(t)  = P[Z(t)  * 1] . 


For  k > 2 , note  that 


(2.58)  k;p[z(t)=ki  5 k:pk(t)  = 


s=0 


Then  (2.58)  applied  to  (2.3)  for  k * 2 yields 


(2.59)  P2(t)  - -j  J h"(l-Q(t-u))P2(t-u)dG(u)  + J*  h'(l-Q(t-u))P2(t 


-u)dG(u) . 
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Since  the  theorem  is  true  for  P(t),  then  as  (2.59)  is  of  the  same 
form  as  (2.8),  a similar  analysis  as  that  used  to  establish  the  result  for 
P(t)  also  establishes  it  for  P2(t).  Assume  the  result  true  for  P^(t),  k < n. 
Applying  (2.58)  to  (2.3)  for  k = nfl  yields  an  equation  of  the  form 

(2.60)  P^Ct)  = ^(t)  + j\'(l-Q(t-u))Pnfl(t-u)dG(u) 
where  the  induction  hypothesis  yields  that 

(2.61)  f(t)  = o(t"4). 

In  view  of  (2.61),  equation  (2.60)  Is  of  the  same  form  as  (2.8),  and 
the  analysis  for  P(t)  applies  to  P .(t),  yielding  the  theorem. 

Remark.  This  result  is  known  for  the  critical  discrete- time  or 
Galton-Watson  process.  See  [2]  C'h.  1,  for  example.  It  may  be  possible 
that  a series  of  comparison  sequences  of  the  type  given  in  (2.44),  (2.45) 
relating  the  critical  Bellman-Harris  process  with  a corresponding  critical 
Galton-Watson  process  could  yield  this  result. 

3.  Total  Progeny. 

For  another  application  of  the  method,  let 

(3.1)  N(t)  =■  total  number  of  progeny  born  by  t in  a critical  age- 
dependent  branching  process  with  •>  h’(l)  sCikPk>0* 

Denote,  for  0 < a < 1,  t > 0 

(3.2)  H(s,t)  * E(sN(t))  - tj^l  P[N(t)«k]a\ 


14 


Then  (e.g.  [4],  p.  394) 


(3.3) 


H(s , t) 


-S 


* 

1 - G(t)  + 


r 

! h(H(s,t-u))dG( 
* 0 


Since,  for  k > 1, 


(3.4) 


, k!P[N(t)=k]  « k!Q  (t) 

I..0 


one  may  apply  (3.4)  to  (3.  3)  to  obtain  the  limiting  values  of 
P[N(t)=k]  s Qk(t). 

For  k = 1,  2 one  obtains,  since  H(0,t)  3 0,  that 


(3.5) 


Qx(t) 


1 - G(t)  + f h(0)dG(u)  - pf 


as  t 


pt/2  Pt 

(3.6)  Q,(t)  - + h ' (0)Q. (t-u)dG(u)  - pQp. 

2 Jt/2  1 u 1 


as  t -*  •. 


Similarly,  by  induction  and  an  application  of  Leibniz'  lemma  for  successive 
differentiation,  the  limiting  values  of  the  fQk(t))  may  be  obtained.  This 
result  is  implied  by  a representation  essentially  as  that  in  [ 1J , pp.  275-276. 
No  assumptions  on  higher  moments  of  G are  required,  but  all  moments  of 
h(s)  are  needed  in  this  approach. 


a 
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